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Computation of the flow field around a projectile moving through a gun 
barrel is a complex problem. The complications are due to presence of 
numerous factors, such as turbulence, compressibility, anisotropy of the 
stress field, complexities due to surface geometry, and change of the 
computational domain in time. In this study, calculation of the flow 
characteristics around 155 mm moving projectile in gun barrel is 
undertaken. The flow is considered to be unsteady, compressible, turbulent, 
and nonisothermal. An anisotropic two-equation model for the Reynolds 
stress that was developed through a novel technique involving the 
representation of the energy spectrum and invariance based scaling is used 
in perfonning calculations. The model was implemented in a commercial 
code (Fluent) and its perfonnance was validated with benchmark 
experimental flows. Computations are performed for stagnation pressures of 
up to 600 MPa, temperatures of up to 5000K, and free-stream velocities of 
about 1000 m/s in this study. The results are analyzed and the predictive 
capability of the model is discussed. 


INTRODUCTION 

Prediction of the flow field around a projectile moving through a gun barrel is 
very complex task. The complexities are due to presence of numerous factors, such as 
turbulence, anisotropy of the stress field, complex surface geometry, change in time of 
computational domain, and compressibility. Turbulent flow field past projectiles are 
known to exhibit a wide range of length and scales, and their accurate representation 
and modelling plays crucial role in the projectile’s design. The current trend in 
simulation of turbulent flows for practical application relies heavily on computations 
based on Reynolds Averaged Navier-Stokes (RANS) equations. However, Reynolds 
stress tensor, embedded in these equations, has to be modelled. Most of these models 
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are constructed using intuitive approximation or physical analogy rather than formal 
manipulation of Navier-Stokes equations [1], [2]. This approach often assumes 
turbulence to be characterized by a single length and time scale while in practice 
turbulence, at high Reynolds numbers, is exited by a wide range of time and length 
scales. In engineering practice the most commonly used Reynolds stress models are 
two-equation models due to their simplicity and robustness. However, to be of general 
validity, such models should incorporate the effects of system rotation, swirl, surface 
curvature, anisotropy and compressibility. Two equation models failed to predict 
rotating turbulence due to the break down of two basic assumptions implicit in these 
models- the Taylor approximation and the universal Kolmogorov spectrum [4]. By 
replacing the universal Kolmogorov energy spectrum with a rotation modified energy 
spectrum the eddy viscosity was refonnulated. Based on the modified spectrum and the 
approach used by Reynolds [5], the equation for dissipation rate was also reformulated. 
An extended Kolmogorov phenomenology for compressible turbulence is used to obtain 
compressible constitutive relation for the Reynolds stresses. The interaction of vortical 
and acoustic modes affects the solenoidal component of the energy spectrum which 
characterizes the vortical turbulence and is responsible for the Kolmogorov energy 
cascade. An eddy viscosity and a constitutive relation for compressible Reynolds stress 
are derived using this compressibility corrected solenoidal energy spectrum. A 
consistent anisotropic two-equation model for the Reynolds-stress is developed with the 
aid of invariance principles and available scaling laws and used in the present 
investigation [6], 


MODEL FORMULATION 


A general constitutive relation for Reynolds stress xq is considered in the following 
form: 
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are nonnalized rate of strain tensors. 

A general form of polynomial representation for bq [7]is given as: 
b = ^G A T A 
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Where, T are the integrity bases, 
J ( "=5’;7 ,2) =5*r-rY ' 


r (3) = S* 2 -^-{s* 2 }/;T ,4) = W* 2 ~ty* 2 }l 

T (5) =W*S * 2 -S*W* 2 ;T (6) = W' 2 S* +S*W * 2 -j{s*W* 2 }l 

t (7) = rv*s*rv* 2 -rv* 2 s*rv*;T (8) =s*rv*s* 2 -s* 2 rv*s* 

T (9) =W* 2 S * 2 + S* 2 W * 2 -j{s* 2 tV* 2 }l;T m =W*S* 2 W * 2 -W* 2 S* 2 W* 

Where, indicates the trace, and I = 5jj. 

In explicit algebraic stress models [8], the unknown scalar functions, G x are obtained 
based on second order closure theories. The basic form of the scalar functions, G x , may 
also determined by the use of scaling laws. This procedure is employed using a 
representation of eddy viscosity that includes the effect of rotation [4]: 
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Here, Q m is the angular velocity of the noninertial frame relative to an inertial frame, a & 
2.25; and the coefficients Ce «0.1, C-j «0.2, [4]. The corresponding length scale, l is 

3 

r=c^-/f+c,w;w; 

If length scales in the linear term and higher order terms in the model are dimensionally 
consistent, it follows that the resulting expressions for the scalar functions, G l can be 
written as: 
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The general form of constitutive relations for three-dimensional flows is restricted to 
third order, the lowest order for characterizing three-dimensional features. 

The Reynolds stress, r,/ can be written as: 
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Coefficient calibration using data of [10] gave C <1, «0.23, C (2) ~ 0.01, C* 4) & 1.2, C^’ 1 « 
0.06,and based on axisymmetric expansion and fully-developed pipe flow, gave: C* 3) 
0.1, C* 5) « 0.1 .Transport equation for K and £ are modeled in the following fonn 
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with i/o=5. 6, p =0.138, C 3 ~ 1.53; tj=(K/s) V( ISjjSj) ; P= -T // iS’, / and «/<; ~ a £ ~ 1.39, are 
derived on a double-expansion procedure and RNG to account for anisotropic 
production that include the effects of rotation [9], [4], Similar approach can be used to 
obtain the anisotropic Reynolds stress tensor and the eddy viscosity [11], and it is 
expressed as 
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v r = c„ £i/ji+c,w;w;(i + a‘m;) (8) 

Where M t is turbulent Mach number, and A* is modeling coefficient determined by 
fixed point analysis of compressible homogenous turbulence. Using, again, the general 
polynomial representation of constitutive relation for Reynolds-stress (2) we obtained 
corrected Ty for compressible flow in similar form as (5) with correction by (1 + A*M?)’ 

i=l, 2, 3. The model developed in this way is of general, tensorially invariant form 
based on scaling laws and accounts for the effects of rotation, swirl, surface curvature, 
and compressibility. It had been earlier validated against experimental results for axially 
rotating pipe, rotating channel, mixing layer, etc. 


SIMULATION AND RESULTS 


In this study, calculation of the flow characteristics around 155 mm moving 
projectile with blow-by in gun barrel is undertaken. The flow is considered to be 
unsteady, compressible, turbulent, and nonisothermal. Simulations were perfonned 
using a commercially available Finite Volume Code (Fluent), with the above formulated 
model incorporated into the code via user define function (UDF). Figure 1. presents 
geometry and meshing of the model (using the Fluent pre-processor Gambit). 



Figure 1. Meshed Model of 155 mm projectile with blow-buy in the gun barrel 

As the projectile moves down the barrel formation of the shock wave is expected to 
occur once projectile accelerated enough. Figure 2 presents the Mach number contour 
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Figure 2. Oblique shock around nose of the projectile 

around nose of the projectile. The turbulence model, as it can be seen in above figure, 
was capable of capturing the oblique shock, as well as reflected waves from the gun 
barrel. The pressure contour at the projectile base is presented by Figure 3. As the 
projectile progressed from starting position where initial pressure was 4287 atm., the 
pressure at base of the projectile regressed due to the movement of the projectile. Blow- 
by section of this set up at the narrowest place has the width of the only 125 microns. 
As the projectile starts moving products of combustions expand and rush into this space 
where jump in the velocity occurs and flow turns supersonic. This section diverges as it 
moves toward the nose and as expected from compressible gas theory it accelerates 
which is clearly visible in the Figure 4. 
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Figure 3. Pressure Distribution behind the Projectile 
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Figure 4. Mach Number Contour for the Blow-by 
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CONLUSIONS AND FUTURE WORK 

Anisotropic Reynolds stress model developed using a phenomenological 
approach was implemented into the commercial code Fluent and used for the simulation 
of 155 mm moving projectile with blow-by down the gun barrel. The model preserves 
the two-equation form and meets the invariance constraints and includes the rotation- 
modified energy spectrum and a generalized constitutive representation that does not 
require the use of local equilibrium assumption. Model coefficients were obtained with 
fixed-point analysis of simple flows and using empirical data. Previous validation has 
shown that model is capable of capturing the effects of rotation, curvature and 
compressibility. In this paper the model was used to predict the dominant features of the 
flow field around the moving projectile. Future work will focus on the use of finer mesh 
to accurately capture the turbulence parameters. A detailed comparison with 
experimental data is also planned for further validation of the model coefficients. 
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